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We present results on time-dependent CP asymmetries in neutral B decays to several CP eigen- 
states. The measurements use a data sample of about 88 million Y(4S) — > BB decays collected 
between 1999 and 2002 with the BABAR detector at the PEP-II asymmetric-energy B Factory 
at SLAC. We study events in which one neutral B meson is fully reconstructed in a final state 
containing a charmonium meson and the other B meson is determined to be either a B or B° 
from its decay products. The amplitude of the CP asymmetry, which in the Standard Model is 
proportional to sin2/3, is derived from the decay-time distributions in such events. We measure 
sin2/3 = 0.741 ± 0.067 (stat) ± 0.034 (syst) and |A| = 0.948 ±0.051 (stat) ± 0.030 (syst). The magni- 
tude of A is consistent with unity, in agreement with the Standard Model expectation of no direct 
CP violation in these modes. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The Standard Model of electroweak interactions de- scribes CP violation in weak interactions as a con- 
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sequence of a complex phase in the three-generation 
Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing ma- 
trix Q . In this framework, measurements of CP asymme- 
tries in the proper-time distribution of neutral B decays 
to charmonium final states provide a direct measurement 
of sin2/3 @, where = arg [ -V cd V£/V tA V t i }. 

Observations of CP violation in B° decays were re- 
ported last year by the BABAR || and Belle [Q collabo- 
rations. The PEP-II collider has since delivered an addi- 
tional 63 fb _1 , thereby approximately tripling the data 
sample near the T(4S) resonance. In this Letter we re- 
port a more precise measurement of sin2/3 using the full 
sample of about 88 million BB decays. The BABAR de- 
tector and the measurement technique are described in 
detail in Refs. || and respectively. Changes in the 
analysis with respect to the published result [|3| include 
processing of all data with a uniform event reconstruc- 
tion, a new flavor-tagging algorithm, and the addition of 
the decay mode B° — > ilcK®. 

We reconstruct a sample of neutral B mesons (Bcp) 
decaying to the final states J/i()K®, ip(2S)K®, xaK®, 
r] c K° s , J/^K*°{K*° -> K° s n°), and J/tpK° L . The J/tf> and 
il>(2S) mesons are reconstructed through their decays to 
e + e~ and the ip(2S) is also reconstructed through 

its decay to J/ipir + n~ . We reconstruct \ci mesons in the 
decay mode J/ipj and rj c mesons in the K^K + tt~ and 
K + K~tt° final states 0. The K® is reconstructed in its 
decay to n + n~ (and to 7r°7r° for the J/ipK® mode). We 
examine each event in the Bcp sample for evidence that 
the recoiling B meson decayed as a B° or B° (flavor tag) . 

The proper-time distribution of B meson decays to a 
CP eigenstate with a B° or B° tag can be expressed in 
terms of a complex parameter A that depends on both 
the B°-B° oscillation amplitude and the amplitudes de- 
scribing B° and B° decays to this final state The 
decay rate f+(f-) when the tagging meson is a B°(B°) is 
given by 



of the final state /. Thus, the time-dependent CP asym- 
metry is 



f±(At) 



-|At|/T B a 
4tro 



2 Im A . , . . . 
1 ± - — sin (Am d At) 



1 + |A| 5 
1 + |A| ; 



cos (Am^At) 



, (1) 



where At = t Tec — ftag is the difference between the proper 
decay times of the reconstructed B meson (B roc ) and 
the tagging B meson (B tag ), Tgo is the B° lifetime, and 
Arrid is the B°-B° oscillation frequency. The sine term in 
Eq. [j] is due to the interference between direct decay and 
decay after flavor change, and the cosine term is due to 
the interference between two or more decay amplitudes 
with different weak and strong phases. CP violation can 
be observed as a difference between the At distributions 
of B°- and B°-tagged events or as an asymmetry with 
respect to At = for either flavor tag. 

In the Standard Model, A = rjfe~ 2113 for charmonium- 
containing b — > ccs decays, where r)f is the CP eigenvalue 



f+(Ai) + f_(Af) 



= — Tjf sin2/?sin (Ato^ At), 



(2) 



with r) f = -1 for J/i>K% ^{2S)K G S , XaK° s , and rj e K°, 
and +1 for J/tpK®. Due to the presence of even (L=0, 
2) and odd (L=l) orbital angular momenta in the B — > 
J/4>K*° final state, there can be CP-even and CP-odd 
contributions to the decay rate. When the angular infor- 
mation in the decay is ignored, the measured CP asym- 
metry in J/ipK*° is reduced by a factor 1 — 2R±, where 
R± is the fraction of the L=l component. We have mea- 
sured R± = (16.0±3.5)% D, which gives r\ f = 0.65±0.07 
after acceptance corrections in the J/ipK*° mode. 

The event selection, lepton and K ± identification, and 
J/tp and ip(2S) reconstruction used in this analysis are 
similar to those described in Ref. |@), as are the selec- 
tion criteria for the channels J/ipKf, ip{2S)Kg, XciK°, 
J/ipK*°, and J/ipK° L . The B° -► 7] C K° S sample selec- 
tion is described in Ref. |l0|. In brief, the K ± candi- 
dates must satisfy kaon identification criteria and the 
— > 7T + 7r _ and ir° — -» 77 candidates are required to 
have reconstructed masses within 12.5 and 15 MeV/c 2 , 
respectively, of their nominal masses JIl| . The r\ c can- 
didates (with 2.90 < Mkktt < 3.15GeV/c 2 ) are com- 
bined with K® — > 7T + 7r _ candidates reconstructed within 
10 MeV/c 2 of the K® nominal mass to form a B° candi- 
date. This sample includes a contribution of (15 ± 2)% 
from hadronic J/ifj decays to the KKir final states. 

We select candidates in the B° -> J/ipK° s , ip{2S)K° s , 
XdK^, and J/tpK*° sample by requiring that the dif- 
ference AE between their energy and the beam energy 
in the center-of-mass frame be less than three stan- 
dard deviations from zero. The AE resolution is about 
10 MeV, except for the mode with K° -> tt°tt° (33 MeV) 
and with K*° (20 MeV). The B° -> rj c K° candidates 
are required to have \AE\ less than 40 (70) MeV for 
the KgK +/ K~ (K + K~7r°) modes. For all modes ex- 
cept J/ifjK®, the beam-energy substituted mass mss = 
V^bcam) - (pT) 2 must be greater than 5.2 GeV/c 2 . 
To determine numbers of events and purities, a signal re- 
gion 5.270 (5.273) < to es < 5.290 (5.288) GeV/c 2 is used 
for modes containing K° s (K*°). In the J/ipK® mode, the 
AE resolution is 3.5 MeV (after B mass constraint) and 
the signal region is defined by \AE\ < 10 MeV. 

A measurement of Aqp requires a determination of 
the experimental At resolution and the fraction w of 
events in which the tag assignment is incorrect. This 
mistag fraction reduces the observed CP asymmetry by 
a factor 1 — 2w. Mistag fractions and At resolution 
functions are determined from a sample of neutral B 
mesons that decay to flavor eigenstates (Bfl av ) consist- 
ing of the channels D^~h + (h + = n + , p + , and ) and 
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J/tp K*°(K*° — ► if+7r ). Validation studies are per- 
formed with a control sample of B + mesons decaying 
to the final states J/ipK^+, ip(2S)K+, XaK + , r/ c K+, 

and D W V. 

We use multivariate algorithms to identify signatures 
of B decays that determine the flavor of -Btag- Primary 
leptons from semileptonic B decays are selected from 
identified electrons and muons as well as isolated ener- 
getic tracks. We use the charges of identified kaon candi- 
dates to define a kaon tag. Soft pions from D* + decays 
are selected on the basis of their momentum and direction 
with respect to the thrust axis of £>tag- A neural network, 
which combines the outputs of these physics-based algo- 
rithms, takes into account correlations between different 
sources of flavor information and provides an estimate of 
the mistag probability for each event. 

By using the outputs of the physics-based algorithms 
and the estimated mistag probability, each event is as- 
signed to one of four hierarchical, mutually exclusive tag- 
ging categories. The Lepton category contains events 
with an identified lepton, and a supporting kaon tag if 
present. Events with a kaon candidate and soft pion with 
opposite charge and similar flight direction are assigned 
to the Kaon I category. Events with only a kaon tag are 
assigned to the Kaon I or Kaon II category depending on 
the estimated mistag probability. The Kaon II category 
also contains the remaining events with a soft pion. All 
other events are assigned to the Inclusive category or 
excluded from further analysis based on the estimated 
mistag probability. The tagging efficiencies e< for the four 
tagging categories are measured from data and summa- 
rized in Table |. The figure of merit for tagging is the 
effective tagging efficiency Q = ^ £^(1 — 2u>i) 2 . This 
algorithm improves Q by about 7% (relative) over the 
algorithm used in Ref. (|]. 

The time interval At between the two B decays is cal- 
culated from the measured separation Az between the 
decay vertices of B lec and B tag along the collision (z) 
axis 0. We determine the z position of the B Tec vertex 
from its charged tracks. The B tag decay vertex is deter- 
mined by fitting tracks not belonging to the B ICC can- 
didate to a common vertex, employing constraints from 
the beam spot location and the B rec momentum |J . We 
accept events with a At uncertainty of less than 2.5 ps 
and |Ai| < 20 ps. The fraction of events satisfying these 
requirements is 95%. The r.m.s. At resolution for 99.7% 
of these events is 1.1 ps. 

The signal region contains 2641 events which satisfy 
the tagging and vertexing requirements. In Table |n] 
we list the number of events and the signal purity for 
the tagged Bqp candidates. The purities are deter- 
mined from fits to the toes (all modes) or AE {K® 
mode) distributions in data, or from Monte Carlo simu- 
lation (K*° mode). Figure [l] shows the toes distribu- 
tion for modes containing a K® or K*° and AE dis- 



TABLE I: Efficiencies e^, average mistag fractions Wi, mistag 
fraction differences Awi = Wi(B°) — Wi(B°), and Q extracted 
for each tagging category i from the _Bfl av and Bcp samples. 



Category 


e (%) 


w (%) 


Aw (%) 


Q (%) 


Lepton 


9.1 ± 0.2 


3.3 ± 0.6 


-1.5 ± 1.1 


7.9 ± 0.3 


Kaon I 


16.7 ± 0.2 


10.0 ± 0.7 


-1.3 ± 1.1 


10.7 ± 0.4 


Kaon II 


19.8 ± 0.3 


20.9 ± 0.8 


-4.4 ± 1.2 


6.7 ± 0.4 


Inclusive 


20.0 ± 0.3 


31.5 ± 0.9 


-2.4 ± 1.3 


2.7 ± 0.3 


All 


65.6 ± 0.5 






28.1 ± 0.7 



tribution for the J/ipK® candidates. For all modes ex- 
cept f] c Kg and J/ipK®, we use simulated events to es- 
timate the fractions of events in the Gaussian compo- 
nent of the toes fits due to cross-feed from other de- 
cay modes. For the r/ c K® mode the cross-feed fraction 
is determined from a fit to the Mkktt and Toes distri- 
butions. These fractions range from (0.3 ± 0.1)% for 
J/ipK° (K° s -> 7T+7T-) to (13.1 ± 5.9)% for i lc K°. For 
the J/ipK® and J/ipK* decay modes, the composition, 
effective 77 j , and AE distribution (J/ipK^ only) of the in- 
dividual background sources are determined either from 
simulation (for B — > J/ip X) or from the to^+^- sidebands 
in data (for fake J/ip -> 

We determine sin2/3 with a simultaneous unbinned 
maximum likelihood fit to the At distributions of the 
tagged Bqp and Ba av samples. In this fit the At distri- 
butions of the Bcp sample are described by Eq. |l| with 
|A| = 1. The At distributions of the -Bfl av sample evolve 
according to the known frequency for flavor oscillation in 
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FIG. 1: Distributions for Bcp candidates satisfying the 
tagging and vertexing requirements: a) rriBS for the final 
states J/tl>K%, ^(2S)K° S , XdK° s , v c K° s , and J/ipK*°(K*° -> 
Kin ), and b) AE for the final state J/ipK ,. 



TABLE II: Number of events iVtag in the signal region af- 
ter tagging and vertexing requirements, signal purity P, and 
results of fitting for CP asymmetries in the Bcp sample and 
in various subsamples, as well as in the -Bfl av and charged B 
control samples. Errors are statistical only. 



Sample 


Atag 


P(%) 


sin2,3 


J^Klrt{2S)Ki, X ciK», nc 

J/tbK^ (nr = + 1) 

MK*°(K*° -> K%°) 

i J J Vy ± \. \ _Z V 7 J. V c / 1 I 


Kg 1506 
988 
147 


94 
55 
81 


0.76 ± 0.07 
0.72 ± 0.16 

22 + 52 


Full CP sample 


2641 


78 


0.74 ± 0.07 


J/^K S , tp(2S)K s , XdK s , 


VcKg only (ri f = -1] 




Till) K" (K U — » -n- + TT~\ 
J/TpK° s (K% -» 7T 7r°) 

XciK° s 


974 
170 
150 
80 


97 
89 
97 
95 


n S9 + n ns 

u.o^ m u.uo 

0.39 ± 0.24 
0.69 ± 0.24 
1.01 ± 0.40 


VcK° s 


132 


73 


0.59 ± 0.32 


Lepton category 
Kaon I category 


220 
400 


98 
93 


0.79 ± 0.11 
0.78 ± 0.12 


Kaon II category 


444 


93 


0.73 ± 0.17 


Inclusive category 


442 


92 


0.45 ± 0.28 


B" tags 


740 


94 


0.76 ± 0.10 


B° tags 


766 


93 


0.75 ± 0.10 


-Bfl av sample 


25375 


85 


0.02 ± 0.02 


B + sample 


22160 


89 


0.02 ± 0.02 



B mesons. The observed amplitudes for the CP asym- 
metry in the Bcp sample and for flavor oscillation in the 
Bfia. v sample are reduced by the same factor 1 — 2w due 
to flavor mistags. Events are assigned signal and back- 
ground probabilities based on the tcies (all modes ex- 
cept J/ipK*° and J/^K°) or AE (J/ipK° L ) distributions. 
The At distributions for the signal are convolved with 
a common resolution function, modeled by the sum of 
three Gaussians Iq]. Backgrounds are incorporated with 
an empirical description of their At spectrum, contain- 
ing prompt and non-prompt components convolved with 
a resolution function || distinct from that of the signal. 

There are 34 free parameters in the fit: sin2/3 (1), 
the average mistag fractions w and the differences Aw 
between B° and B° mistag fractions for each tagging 
category (8), parameters for the signal At resolution 
(8), and parameters for background time dependence 
(6), At resolution (3), and mistag fractions (8). We fix 
t b o = 1.542 ps and Am d = 0.489 ps^ 1 Q. The de- 
termination of the mistag fractions and At resolution 
function parameters for the signal is dominated by the 
high-statistics i?fl av sample. The measured mistag frac- 
tions are listed in Table ffl. Background parameters are 
determined from events with tties < 5.27GeV/c 2 (except 
J/ipK® and J/ipK*°). The largest correlation between 
sin2/3 and any linear combination of the other free pa- 
rameters is 0.13. We observe a bias of 0.014 ± 0.005 in 
the fitted value of sin2/3 in simulated events. Part of this 
bias (0.004) is due to a correlation between the mistag 
fractions and the At resolution not explicitly incorpo- 



G 

rated in the fit. Therefore we subtract 0.014 from the 
fitted value of sin2/3 in data and include 0.010 in the 
systematic error. 

The fit to the Bcp and Bflav samples yields 

sin2/3 = 0.741 ± 0.067 (stat) ± 0.034 (syst). 

Figure ^| shows the At distributions and asymmetries in 
yields between B° tags and B° tags for the rjf = — 1 and 
rjf = +1 samples as a function of At, overlaid with the 
projection of the likelihood fit result. 




-0.5 - 



-5 5 

At (ps) 

FIG. 2: a) Number of rjf = -I candidates (J/ipKs, 
ip(2S)Kg, XciKg, and rj c Ks) in the signal region with a B° 
tag N B o and with a B tag Ng , and b) the raw asymme- 
try (N B o - N S o)/(N B o + N S o) as functions of At. The solid 
(dashed) curves represent the fit projection in At for B° (B°) 
tags. The shaded regions represent the background contribu- 
tions. Figures c) and d) contain the corresponding informa- 
tion for the rjf = +1 mode J/ipK^. 



The dominant sources of systematic error are the un- 
certainties in the level, composition, and CP asymme- 
try of the background in the selected CP events (0.023), 
the assumed parameterization of the At resolution func- 
tion (0.017), due in part to residual uncertainties in the 
internal alignment of the vertex detector, and possible 
differences between the -Bfl av and Bcp mistag fractions 
(0.012). The total systematic error is 0.034. Most sys- 
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tematic errors are determined with data and will continue 
to decrease with additional statistics. 

The large Bqp sample allows a number of consistency 
checks, including separation of the data by decay mode, 
tagging category, and B tag flavor. The results of fits to 
these rjf = — 1 subsamples are shown in Table [ij and 
found to be statistically consistent. The results of fits to 
the control samples of non-CP decay modes indicate no 
statistically significant asymmetry. 

We also measure the parameter |A| in Eq. [I] from 
a fit to the r/f = — f sample, which has high purity 
and requires minimal assumptions on the effect of back- 
grounds. This parameter is sensitive to the difference 
in the number of B°- and B°-tagged events. In or- 
der to account for differences in reconstruction and tag- 
ging efficiencies for B° and B° mesons, we incorporate 
five additional free parameters in this fit. We obtain 
|A| = 0.948 ± 0.051 (stat) ± 0.030 (syst). The coeffi- 
cient of the s'm(AmdAt) term in Eq. [I] is measured to 
be 0.759 ± 0.074 (stat). The dominant contribution to 
the systematic error for |A|, conservatively estimated to 
be 0.025, is due to interference between the suppressed 
b — > ucd amplitude with the favored b — > cud ampli- 
tude for some tag-side B decays. The other sources of 
systematic error for |A| are the same as in the sin2/3 mea- 
surement. 

This measurement of sin2/3 supersedes our previous re- 
sult Q and improves upon the precision of each of the 
previous measurements || Q by a factor of two. While 
the measured value is consistent with the range implied 
by the measurements and theoretical estimates of the 
magnitudes of CKM matrix elements in the context of 
the Standard Model, it provides a precise and model- 
independent constraint on the position of the apex of the 
Unitarity Triangle 
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